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Abstract. We propose the entropy density as the thermodynamic condition driving 
best the chemical freeze-out in heavy-ion collisions. Taking its value from lattice 
calculations, we find that it is excellent in reproducing the experimentally estimated 
freeze-out parameters. The two characteristic endpoints in the freeze-out diagram are 
reproduced as well. 
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1. Introduction 

At critical temperature T c , the hadronic matter is conjectured to dissolve into quark- 
gluon plasma (QGP). Reducing the QGP temperature disposes hadronization and the 
interacting system goes into chemical equilibrium and finally the produced particles 
freeze out. Below T c h, thermal equilibrium takes place and the compositions of matter 
turn to confine into hadrons. The best way to determine the freeze-out parameters is 
to combine various particle ratios in order to obtain a window in the T c h — \ib diagram 
compatible with the experimental results. \xb is the baryo-chemical potential. The 
question we intent to answer is: what is the universal condition describing the freeze- 
out parameters? [TJ |2| 

Without energy input the chemical reactions always proceed toward chemical 
equilibrium, i.e. balancing particle-absorption and -production. To study the 
consequences of this equilibrium, we recall the equilibrium constant C c h- In a chemical 
reaction like aA + bB — > cC + dD, C c h can be calculated according to the "law of mass 
action" glU] 

C ch = ([C] c ch [D] d ch )/([Ar ch [B] b ch ), (1) 

where (A, B) and (C, D) refer to the reactants and products, respectively, (a, b) and 
(c, d) are the corresponding concentrations. In heavy-ion collisions, the decay channels 
are - in structure - similar to the above chemical reaction. The backward direction can 
be viewed as annihilation/absorption processes. Analogy to Eq. ((JJ we deal so far with 
one specific type of hadron interactions. The hadron system and its approach 
towards the chemical equilibrium are allowed to have many dynamical processes. The 
most important ones are particle-absorption and -production. The different mechanisms 
that drive the system towards the chemical equilibrium are to be taken into consideration 
by summing over all hadron resonances. 

The total free energy of the system SG can be used to determine the likely direction. 
When the reactant and product concentrations are given, then for an ideal gas (no 
enthalpy change), 5G = 5G° + Tln(C c ^). At equilibrium free energy gets minimum (5G 
vanishes) and pressure is constant. Then C c h = exp(— <5G°/T). SG° = 5E — T5S is the 
difference of free energies of products and reactants at standard state, i.e. at T ^ T^. 

C ch ~ exp(-<5s - fi8n/T), (2) 

At equilibrium, the entropy gets maximum. The change in particle number is minimum. 
Therefore, the equilibrium constant is mainly influenced by entropy change. 

2. Thermodynamic conditions 

In the following, we list some thermodynamic expressions for one particle and its anti- 
particle in Boltzmann limit. 

n(T^ B ) = ±Tm*K*(£) sinh(^), (3) 
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where g is the spin-isospin degeneracy factor and Ki are the modified Bessel functions. 
They will be summed over all resonances taken into account. These quantities are related 
to the consequences of chemical equilibrium via Eq. (j2J. In this letter, full quantum 
statistics has been properly taken into account. 

\nZ(T,n B ) = ^ [°°k 2 dk\n\l± e^l, (6) 

27T Z JO L 

where e = (k 2 + m 2 ) 1 ^ 2 is the single-particle energy and ± stand for bosons and fermions, 
respectively. We apply the hadron resonance gas model (HRGM) El 13 E] in order 
to study the conditions driving the freeze-out. All observed resonances up to mass 
2 GeV are included. The particle decays are entirely left away. We use grand canonical 
ensemble and full quantum statistics. Corrections due to van der Waals repulsive 
interactions and the excluded volume have not been taken into account. 
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Fig. 1. Freeze-out curve under the condition of constant energy per particle The 
points are the parameters taken from indicated accelerators. 



Before we report the results, we remind of the two characteristic points of the freeze- 
out diagram; One at T ch ^ and fiB = and the other at T ch = and very large \ib- 
Localizing the first point has been the subject of different experimental studies 
It has been found that T c h(f^B = 0) ~ 174 MeV ^2]- The lattice estimation for the 
deconfinement temperature is T c (fiB = 0) = 173 ± 8 MeV ^3j- This implies that the 
deconfinement and freeze-out lines seem to be coincident at low /xg. For the second 
point, we are left with applying effective models. As T — > 0, the nucleons in the hadron 
gas get dominant. Applying Fermi statistics, Eq. EJ we find that the chemical potential 
corresponding to the normal nuclear density, no ~ 0.17 fm -3 is fi c h ~ 0.979 GeV. This 
value can slightly be different according to the initial conditions |14j . 
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Fig. 2. Freeze-out curve under the condition ub = 0.17 fm~ 3 9 . In calculating the 
short-dashed curve |10| . the percolation theory has been taken into account. 



A proposal to describe the freeze-out parameters at different incident energies has 
been reported in Ref. j^j. The authors started from phenomenological observations 
at GS1/SIS energy that e/n k 1 GeV. They have applied Boltzmann approximations 
in calculating e/n, Eq. (JEJ), and analytically handled the hadron resonance gas at 
T ps 50MeV ~ 0.8 GeV), i.e. the freeze-out parameters at GSI/SIS energy as a 
Fermi gas of degenerate nucleons, Eq. At high T and small the authors assumed 
that the pions and rho- mesons get dominant. The baryons are entirely neglected. 
Applying Eq. (jHJ) with the effective mass m = m n + m p they got almost the same 
value for the ratio e/n 0. In Fig. ^ we use HRGM in order to map out T — fi c h freeze- 
out diagram according to this condition. In calculating n and e, we take into account 
full quantum statistics and all observed resonances up to mass 2 GeV. 

In Ref. jHj, total baryon number density % was imposed to interpret the chemical 
freeze-out curve. For this scenario, the baryon-baryon and baryon-meson interactions 
were assumed to drive the chemical equilibrium. For the value n% = 0.12 fm -3 , two- 
third the normal nuclear density, it has been argued that it depends on the correction 
due to van der Waals repulsive interactions. In Fig. El we plot HRGM results under 
the condition ub = 0.17 fm~ 3 , the normal nuclear density We use this value, since we 
entirely leave away all corrections. We find that our results suggest that the assumption 
of baryon-baryon and baryon-meson and the repulsive interactions are not well-founded. 
Nevertheless, we find that the value ns = 0.17 fm -3 is satisfactorily able to reproduce 
the two endpoints of freeze-out diagram. 

Another model we consider is in the framework of percolation theory jTOj. At T = 0, 
the freeze-out occurs when the nucleons no longer form interconnected matter. The 
corresponding density is found to be ~ 0.17 fm -3 and consequently, fiB = 0.979 GeV. 
At a B = 0, it has been found that T c h ~ 175 MeV for 7 S = 0.5 [10J. 7 S gives the 
strangeness saturation. The results are shown in Fig. El 

We can so far conclude that the last two models |H1 CHj are able to reproduce the 
two endpoints. Both apparently overestimate the freeze-out parameters at BNL/AGS 
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and GSI/SIS energies. Model |2] describes well the freeze-out parameters at low energy. 
It slightly underestimates the parameters at BNL/RHIC and CERN/SPS energies. Its 
largest discrepancy is at energies lower than GSI/SIS. At \lb = 0.979 GeV, which 
as given above corresponds to n , we find that the freeze-out under the condition 
e/n = 1 GeV occurs at T c h ~ 35 MeV! Furthermore, extrapolating the resulting curve 
to abscissa results in particle number density 25 — 30 n . According to Eq. [7[ the ratio 
e/n at T = and \l b = 0.979 GeV equals 2.89 GeV. 
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This relation is valid for fermions and therefore, we drop out the exponent ±1 [To] . 
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Fig. 3. Lattice QCD results on entropy density normalized to T 3 for 2 and 2 + 1 
quark flavors at \±b = 0. The thin curves give HRGM results for physical masses. The 
results for re-scaled resonance masses are given by the other two curves. 
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Fig. 4. The freeze-out curve under the condition of constant s/T 3 . For non-strange 
resonances, s/T 3 = 5 and for all resonances, s/T 3 = 7 (Fig.[3J). The two solid circles 
as well as the experimentally estimated points are very well reproduced. 
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Fig. 5. T c h, s and e/n are calculated along the freeze-out curve. We notice that s 
decreases much faster than T c h, so that at very large ug the condition s/T 3 = 7 is no 
longer valid. The ratio e/n. = 1 GeV does not remain constant along /i^-axis. 

3. Entropy density at chemical freeze-out 

As mentioned above, our condition for chemical equilibrium in heavy-ion collisions is 
the entropy density, Eq. 4\ . For vanishing free energy, the equilibrium entropy gives 
the amount of energy which can't be used to produce additional work. In this context, 
the entropy can be defined as the degree of sharing and spreading the energy inside the 
system. The way of distributing the energy is not just an average value, but rather the 
way that controls the chemical equilibrium. 

As T — > and /ig ^ 0, the thermodynamic quantities are given in Eq. (JJJ). In this 
limit, the entropy density vanishes and HRGM is no longer applicable. At //# = 0, the 
system becomes meson-dominant. In this case, the entropy is finite [HJG31- What is 
the entropy that characterizes the freeze-out at finite T and yUe? As mentioned above, 
the deconfinement transition is coincident with the freeze-out at small \xb- We therefore 
rely on the lattice calculations El El CHI- in Fig- El we plot the lattice results on 
s/T 3 vs. T/T c at /i£ = [El E] for rif = 2 and rif = 2 + 1 quark flavors. To compare 
with the lattice calculations in which very heavy quark masses are used, the resonance 
masses have to be re-scaled to values heavier than the physical ones [5J E] . The results 
with the physical masses are given by the thin curves. The two horizontal lines point at 
the value of s/T 3 at corresponding T c . For physical masses, we find that s/T 3 = 5 for 
rif = 2 and s/T 3 = 7 for rif — 2 + 1. The normalization with respect to T 3 should not 
be connected with massless ideal gas. Either the resonances in HRGM or the quarks on 
lattice are massive. 

At given we calculate T c h according to constant s/T 3 . The results are plotted 
in Fig. |U The dotted curve represents rif = 2 results, (s/T 3 = 5). The solid curve gives 
rif = 2 + 1 results, (s/T 3 = 7). We find that the characteristic endpoints as well as all 
experimentally estimated freeze-out parameters are very well reproduced. Comparing 
with Fig. and Fig. |2j one finds that our results fit best freeze-out parameters. In 
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Fig. |5J we plot e/n, T ch and s in dependence on \i B along the freeze-out curve. We 
notice that as \1b increases, both s and T ch decrease, too. We also find that s decreases 
much faster than T. The ratio s/T 3 becomes greater than 7 at very large \xb- In this 
limit, HRGM is no longer applicable and the numerics terminates. At very large 
thermal entropy s is expected to vanish, since it becomes proportional to T (third law of 
thermodynamics). The quantum entropy [T7[ [TBI EE 1201 I2H 122] is entirely disregarded 
here. We also find that e/n does not seem to remain constant along the freeze-out curve 
as assumed in [2.. Since HRGM is not applicable in the limit T = 0, we apply Fermi 
statistics, Eq. as given above. 

4. Conclusion 

We reviewed the conditions suggested to describe the freeze-out parameters. We used 
HRGM in order to map out the freeze-out curve according to these conditions. Full 
quantum statistics has been properly taken into account. We compared the results of the 
models proposed and check their abilities in reproducing the experimentally estimated 
freeze-out parameters and the characteristic endpoints (Fig.QJand Fig.|2J). We found that 
there are different constraints in reproducing the endpoints and fitting the parameters. 

We proposed the entropy density s as the thermodynamic condition driving the 
chemical freeze-out, Eq. |21 Taking its value from lattice QCD simulations at [is — 
and assuming that s normalized to T 3 remains constant in the entire /x^-axis, we 
obtained the results shown in Fig. HJ The freeze-out parameters T ch and \ib are 
very well described under this condition. The two characteristic endpoints are also 
reproduced. We conclude that the given ratio s/T 3 characterizes very well the final states 
observed in all heavy-ion experiments. Increasing the incident energy leads to increasing 
particle yields. The production rates of particles decrease exponentially with their 
masses. This phenomenological observations are regarded in the way, that increasing 
energy/temperature is considered by including heavy resonances. Changing energy 
with changing particle number is given by chemical potential \xb- The amount of 
energy which produces no additional work, i.e. vanishing free energy, is the entropy at 
chemical equilibrium. 

I thank Rudolf Baier, Tamas Biro, Rajiv Gavai, Krzysztof Redlich and Helmut Satz 
for useful discussions. 
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